Implications of reflectance measurements on the mechanism for superconductivity in 

MgB 2 
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Recent optical studies in c-axis oriented superconducting MgB2 films indicate that the electron- 
phonon coupling is weak 111 . We reinforce this conclusion by examining the raw reflectance data; 
its frequency dependence is incompatible with strong electron-phonon scattering. This is further 
strengthened by analysis of the real part of the conductivity, and by the temperature dependence 
of the effective Drude scattering rate. Using a realistic electron-phonon spectral shape J2|, we find 
Ati ~ 0.15, in agreement with Tu et al. M. To the extent that Atr ~ A, this disagrees sharply with 
model calculations H-ffl, and is far too low to provide the means for T c = 39 K. A simple model 
is constructed with coupling to a high frequency excitation, which is consistent with both the low 
frequency optical data and the high T c . 



The discovery of superconductivity in MgB 2 at 39 K 
H has naturally prompted questions about mechanism. 
Thus far calculations of the electron-phonon coupling 
strength 0-B] suggest that the conventional electron- 
phonon mechanism may be responsible for T c = 39 K. 
Nonetheless, there remains some concern that a signifi- 
cant modification of the theory is required to explain the 
high T c , due to anharmonic || or non-adiabatic jTjJg] ef- 
fects. In addition, non-electron-phonon mechanisms have 
been proposed P,[l0[. 

All of these deliberations have been essentially theo- 
retical in nature. On the experimental side, the isotope 
effect measurements pd|Jl^ ] can be understood within 
the electron-phonon framework, but only with very high 
frequency phonons, and a very large Coulomb repulsion 
|l2| , [l3| . Very recently reflectance measurements on c-axis 
oriented thin films have been analyzed within the con- 
ventional electron-phonon framework W . These authors 
conclude that the electron-phonon coupling in MgP>2 is 
very weak, X tr = 0.13 ± 0.02. The purpose of this paper 
is to further analyze their data, with the goal of explor- 
ing possible ways in which the reflectance data can be 
interpreted with a more sizeable electron-phonon inter- 
action. We find, in agreement with Tu et al. jjj that the 
reflectance data are not consistent with an electron- 
phonon interaction sufficiently strong to produce T c = 39 
K. Consistency can be achieved if another excitation is 
assumed to exist at much higher frequency. 

Through Kramers-Kronig (KK) analysis of the re- 
flectance data, Tu et al. obtain a very Drude-like 
response at low frequencies. As noted previously for 

Bai^K^BiOs (BKBO) @-@ this wil1 be the case even 
if the electron-phonon interaction is sizeable. In fact, 
the low frequency Drude result for the real part of the 
conductivity is given by (see Eq. (13) in p"5|) 
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where, at low temperature, l/r c ff = l/ri mp /(l + A tr ), 
where l/r; m p is the elastic scattering rate due to impu- 
rities and Atr is the mass enhancement parameter de- 
fined specifically for transport p^H. Similarly wp G fF = 
ujp/ \/l + Atr is the effective plasma frequency, with wp 
the bare plasma frequency. This expression properly re- 
duces to the correct zero frequency limit, where renormal- 
ization effects are known to cancel fl9[| . In what follows 
we make no distinction between the conventional mass 
enhancement factor A and the transport version, At r , al- 
though the latter is to be used in the conductivity while 
the former enters the single electron self energy |l7| , [l8| . 
We will return to this point later. The experimentally 
derived parameters are determined in the normal state 
at T = 45 K 0: w Poff = 13600 cm- 1 , and l/r off = 75 
cm . Any model with higher energy excitations that 
interact with the electrons (such as phonons) must cor- 
rectly reproduce this low frequency behaviour. 

To attempt to model this data within the conventional 
electron-phonon framework we compute the normal state 
self-energy [g(J an d compute the optical conductivity ef- 
fectively in the bubble approximation Jiq , pl| , since ver- 
tex corrections are implicitly contained in the use of At r 
rather than A. Calculation of the complex conductivity 
allows one to compute all the optical functions, and in 
particular the reflectance ]22| . 

In Fig. 1 we show the reflectance as a function of fre- 
quency, for a variety of model calculations with different 
overall electron-phonon coupling strengths. For this pur- 
pose we have used the model electron-phonon spectrum 
of Kong et al. [g| for the spectral function, a 2 F(tu). As 
calculated, it has A = 0.85, but we have scaled it by a 
constant to achieve any desired coupling strength. It is 
characterized by one dominant peak at u = 74 meV (600 
cm -1 ), which manifests itself in the reflectance as a rel- 
atively sharp downward trend, the magnitude of which 
depends on the coupling strength. Even without the KK 
analysis, then, one can see directly from the reflectance 



data that the electron-phonon coupling in MgB2 must 
be weak. This statement is independent of the precise 
form of the spectral function used, as we have verified 
with other model spectra. The key point is that over the 
frequency range of the phonons (and a little beyond) the 
frequency-dependent scattering rate increases in propor- 

and this results in 



tion to the coupling strength 
a reduced reflectance in this frequency range. In this re- 
spect we should mention that the model calculation with 
A = 0.15 fits the data shown in Fig. 1 very well, but 
nonetheless fails at higher frequency. Specifically the re- 
flectance data [Q is significantly lower than the computed 
result, and shows that at higher frequencies other aspects 
of the problem are required [gij . 

Nonetheless, a model which both fits the low frequency 
reflectance and is capable of producing T c — 39 K con- 
sists of a weak electron-phonon spectrum, along with a 
high frequency excitation (of unknown origin) [ |25| , p6[ , as 
shown in Fig. 2. For simplicity we have adopted a broad 
Lorentzian spectrum centred at 0.5 eV. The coupling 
strength to this additional peak is weak (A « 0.4), but, 
because it is at high frequency, the spectrum shown in 
Fig. 2 results in T c = 39 K (with fi*(uj c = beV) = 0.17). 
With this spectrum the optical reflectance is calculated 
and shown as the dotted curve in Fig. 1; it is barely 
distinguishable from the previous result with A = 0.15. 
Note, however, that we have used A = 0.2 for the scaled 
Kong spectrum, because the high frequency Lorentzian 
contribution essentially renormalizes the low frequency 
part: a 2 ca F{u) ~ a 2 F(uj)/(l + A Lor ) |||. Such a spec- 
trum also improves the agreement with the higher fre- 
quency reflectance data, but we do not attach a great deal 
of significance to this fact since other factors also con- 
tribute |24|. We therefore have not made any attempts 
to optimize this agreement. 

Note that the calculations shown in Fig. 1 all utilize 
different values of the plasma frequency and the impurity 
scattering rate. For example, for the curve calculated 
with A tr = 0.85, we used lu p = 13600-y/l + A tr = 18470 
cm -1 and 1/T imp = 75(1 + A tr ) = 138 cm -1 . 

In Fig. 3 we illustrate the same point through the real 
part of the optical conductivity. In Fig. (3a) we show 
the various models already discussed in association with 
Fig. 1. In Fig. (3b) we focus on the 'Holstein' sideband 
near and above the characteristic phonon energies. It is 
clear that a strong electron-phonon coupling (i.e. of order 
unity or more) disagrees with the data. Moreover, a two- 
band model would not immediately remedy this problem 
since this would lead to extra absorption, which would 
contradict the data. Also shown (dot-dashed curve) is 
the calculation with the spectrum of Fig. 2. Because of 
the renormalization of the low frequency spectrum dis- 
cussed above, the result agrees very closely with the weak 
coupling calculation. 

Another way in which one can analyze the data is 
through the temperature dependence of some of the op- 
tical functions. For this purpose we focus on the tem- 
perature dependent scattering rate P3,f5[ , as defined by 



Eq. (14) of the latter reference. This scattering rate was 
used to fit the low frequency optical conductivity with 
a Drude model at various temperatures JJJ. Two points 
determined experimentally are indicated in Fig. 4. In 
addition we have plotted the results expected for vari- 
ous electron-phonon coupling strengths, again using the 
Kong et al. Q spectrum as a model for the spectral func- 
tion. The data indicates that the effective scattering rate 
more than doubles between 45 K and 295 K; this shows 
that some inelastic scattering must be present. However, 
the required coupling strength to reproduce this temper- 
ature dependence is At r ~ 0.15. This is in agreement with 
our analysis of the frequency-dependent reflectance, and 
agrees with the estimates of Tu et al. based on the optical 
sum rule and the temperature dependent DC resistivity. 
As they pointed out, this coupling strength is much too 
small to account for electron-phonon driven superconduc- 
tivity in this material. Our ad hoc model of Fig. 2 again 
reproduces the correct behaviour (dot-dashed curve), in- 
dicating that some unknown high frequency excitation 
may be responsible for superconductivity. 

In summary we have further analyzed the frequency 
and temperature dependence of the optical reflectance 
and conductivity in c-axis oriented superconducting 
MgB2 films. We have found very strong support for 
the conclusion of Ref . W - — namely the electron-phonon 
coupling strength as derived from optical measurements 
is much too small to account for superconductivity in 
MgB2. In all, five different (but not independent) anal- 
yses have produced this conclusion: 1) the discrepancy 
between the plasma frequency as determined by the low 
frequency Drude fit to the conductivity and that deter- 
mined by the sum rule |3|, 2) the fit to the temperature 
dependence of the DC resistivity data Q, 3) the low fre- 
quency dependence of the reflectance, 4) the magnitude 
of the Holstein side band in the real part of the conduc- 
tivity, and 5) the temperature dependence of the Drude 
width of the low frequency conductivity. These conclu- 
sions have been obtained through calculation with a re- 
alistic electron-phonon spectral shape, obtained through 
calculation gj. Nonetheless, they are quite general, and 
depend only on the fact that the phonons in this mate- 
rial extend up to a maximum of 100 meV. The data of 
Tu et al. H shows very convincingly that the coupling to 
excitations up to these energies must be weak. 

In our opinion the most compelling test is the fre- 
quency dependence of the reflectance, since this prop- 
erty emerges essentially as raw data, and does not re- 
quire any KK analysis. If we accept the possibility of 
a high frequency excitation coupling to electrons to pro- 
duce superconductivity, as is modeled in Fig. 2, then 
other consequences would immediately follow. For ex- 
ample all superconducting properties ought to be very 
BCS-like. 

However, a possiblility for reconciling these results 
with electron-phonon driven superconductivity is that At r 
and A differ significantly. This would be the case, if, for 
example, the electron-phonon scattering was dominated 



by forward scattering |2£| . In this case one should expect 
other unusual normal and superconducting state prop- 
erties; this possibility should be further explored, both 
experimentally and theoretically. 



ACKNOWLEDGMENTS 

I would like to thank Jorge Hirsch for bringing these 
experimental results to my attention, and for subsequent 
discussions. This work was supported by the Natural 
Sciences and Engineering Research Council (NSERC) of 
Canada and the Canadian Institute for Advanced Re- 
search. 



[7] 
[8] 



[12; 

[13 

[14 



[15 

[16 

[17 



J.J. Tu, G.L. Carr, V. Perebeinos, C.C. Homes, M. Stron- 
gin, P.B Allen, W.N. Rang , E.M. Choi, H.J. Kim, and 
Sung-Ik, |cond-mat/0107349| (2001). 

Y. Kong, O.V. Dolgov, O. Jepsen, and O.K. Andersen , 
Phys. Rev. B 64, 020501 (R) (2001); |cond-mat/0102499 



J. Kortus et al, Phys. Rev. Lett. 86 4656 (2001); cond 



mat/0101446 (2001) 



J.M. A n and W.E. Picket t, Phys. Rev. Lett. 86 4366, 
(2001); 



:ond-mat/0102391 (2001) 



J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, 

and J. Akimitsu, Nature 410, 63 (2001). 

T. Yildirim, O. Gulseren, J. W. Lynn, C. M. Brown, T. 

J. Udovic, H. Z. Qing, N. Rogado, K.A. Regan, MA. 

Hayward, J.S. Slusky, T. He, M.K. Haas, P. Khalifah, 

K. Inu maru, and R.J. Cava , Phys. Rev. Lett. 87, 037001 

(2001); 



ond-mat/0103460 (2001) 



A.S. Alexandrov, cond-mat/0104413 (2001). 



E. Cappellu ti, S. Ciuchi, C. Gr imaldi, L. Pietronero, and 
S. Strassler hond-mat/010535l| (2001). 



J.E. Hirsch, Phys. Lett. A282, p. 392 (2001), cond 



mat/0102115 



(2001); J.E. Hirsch and F. Marsiglio, |cond- 
(2001) 



mat/0102479| . 

M. Imada, |cond-mat/0103006| (2001) 



S.L. Bud'ko, G Lapertot, C. Petrovic, C.E. Cunning- 
ham, N. Anderson, and PC. Canfield, Phys. Rev. Lett. 
86 1877 (2001). 

D.G Hinks, H. Claus, and J.D . Jorgensen, Nature 411 
457 (2001) dcond-mat/0104242L. 



A. Knigavko and F. Marsiglio, |cond-mat/0105512j (2001). 
O.V. Dolgov, E.G. Maksimov and S.V. Shulga, in 
Electron- Phonon Interaction in Oxide Superconductors, 
edited by R. Baquero (World Scientific, Singapore 
(1991)),p. 30; S.V. Shulga, O.V. Dolgov and E.G. Mak- 
simov, Physica C 178, 266 (1991). 

F. Marsiglio and J.P Carbotte, Phys. Rev. B 52, 16192 
(1995). 

F. Marsiglio, J.P. Carbotte, A. Puchkov and T. Timusk, 
Phys. Rev. B53, 9433 (1996). 
H. Scher, Phys. Rev. Lett. 25, 759 (1970). 



[18] P.B. Allen, Phys. Rev. B3 305 (1971). 

[19] R.E. Prange and L.P. Kadanoff, Phys. Rev. 134 A 566 
(1964). 

[20] P.B. Allen and B. Mitrovic, in Solid State Physics, edited 
by H. Ehrenreich, F. Seitz, and D. Turnbull (Academic, 
New York, 1982) Vol. 37, p.l. 

[21] F. Marsiglio and J.P. Carbotte, |cond-mat/0106135 
(2001). 

[22] T. Timusk and D.B. Tanner in Physical Properties of 
High Temperature Superconductors I, edited by D.M. 
Ginsberg (World Scientific, Singapore, 1989) p. 339. D.B. 
Tanner and T. Timusk in Physical Properties of High 
Temperature Superconductors III, edited by D.M. Gins- 
berg (World Scientific, Singapore, 1992) p. 363. 
F. Marsiglio and J.P. Carbotte, Aust. J. Phys. 50, 975 
(1997); Aust, J. Phys. 50, 1011 (1997). 
For example, it is clear from other data that other bands 
are relevant at these energy scales, and anisotropic scat- 
tering also complicates the analysis. 

F. Marsiglio and J.P. Carbotte, Phys. Rev. B 39, 2726 
(1989). 

H.J. Kaufmann, O.V. Dolgov, and E.K.H. Salje, Phys. 
Rev. B58 9479 (1998). 

A.V. Puchkov, T. Timusk, W.D. Mosley and R.N. Shel- 
ton, Phys. Rev. B50 4144 (1994). 
ML. Kulic, Physics Reports 338, 1 (2001). 



1.00 



[23] 
[24] 

[25] 
[26] 
[27] 
[28] 



x 



0.95 



0.90 



0.85 



0.80 




A. = 0.15 
I = 0.85 
1 = 2.0 

A, = 0.2 + Lor 
data * 



u eff 



13 600 cm 



1/x eff = 75cm 



-1 







500 



1000 
<B(crrf 



1500 



2000 



FIG. 1. Reflectance vs. frequency; experimental results are 
shown with symbols, and curves indicate calculations, tak- 
ing into account different amounts of electron-phonon scat- 
tering. The plasma frequency and the impurity scattering 
rate are determined for each coupling strength to yield the 
effective plasma frequency and the effective scattering rate 
indicated. The data is in clear agreement with very weak 
electron-phonon coupling, including the case where a higher 
frequency excitation also couples to the electrons. 
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FIG. 2. The model spectral function, consisting of two 
pieces: the electron-phonon part, taken from Kong et al. M 
and scaled to give A = 0.2, and a Lorentzian shaped contri- 
bution peaked at 0.5 eV, and of unknown origin. 




FIG. 4. The temperature dependent effective scattering 
rate, determined by the amount of electron-phonon coupling. 
The two data points reported are indicated with symbols. 
They also indicate that the electron-phonon coupling is weak. 
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FIG. 3. Real part of the optical conductivity vs. frequency, 
(a) showing the Drude-like peak that has been measured in 
thin film samples, and (b) on a finer scale, emphasizing the 
Holstein absorption due to the electron-phonon coupling. As 
in Fig. 1, the plasma frequency and impurity scattering rate 
were determined separately for each coupling strength to yield 
the observed low frequency Drude-like peak. The data (not 
shown) are in much better agreement with the weak elec- 
tron-phonon coupling scenario. 



